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Abstract. Theoretically, the phase transition between the confined and deconfined 
phases of quarks can have a remarkable effect on the spin properties of millisecond 
pulsars and on the spin distribution of the population of x-ray neutron stars in low- 
mass binaries. In the latter class of stars, the effect has already been observed — a 
strong clustering in the population in a narrow band of spins. The observed clustering 
cannot presently be uniquely assigned to the phase transition as cause. However, there 
is another possible signal — not so far observed — in millisecond pulsars that we also 
discuss which would have the same origin, and whose discovery would tend to confirm 
the interpretation in terms of a phase transition in the stellar core. 



1 Motivation 

One of the great fascinations of neutron stars is the deep interior where the 
density is a few times larger than the density of normal nuclei. There, in matter 
inaccessible to us other than fleetingly in relativistic collisions, unfamiliar states 
may exist. The most exotic of these is, of course, quark matter — the deconfined 
phase of hadronic matter. It is quite plausible that ordinary canonical pulsars — 
those like the Crab, and more slowly rotating ones — have a quark matter core, or 
at least a mixed phase of quark and confined hadronic matter. If one were able 
to close pack nucleons to a distance such that they were touching at a radius 
of 0.5 fm, the density would be a mere l/[(2r)^/9o] = 6.5 larger than normal 
nuclear density. If nucleons were close packed to their rms charge radius of 0.8 
fm, the density would be 1.6 times nuclear density. Of course, fermions cannot be 
close packed according to the exclusion principle since their localization would 
give them enormous uncertainty in their momentum. They would be torn apart 
before they could be squeezed so much. This simple, and perhaps oversimplified 
argument, makes it plausible that ordinary (slowly rotating pulsars) have a quark 
matter core essentially from birth. 

By comparison, millisecond (ms) pulsars are centrifugally flattened in the 
equatorial plane and the density is diluted in the interior. We shall suppose that 
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the critical phase transition density lies between the diluted density of millisec- 
ond pulsars and the high density at the center of canonical pulsars. Then as a 
millisecond neutron star spins down because it is radiating angular momentum 
in a broad band of electromagnetic frequencies as well as in a wind of particle- 
antiparticle pairs, or as a canonical neutron star at some stage begins accreting 
matter from a companion and is spun up, a change in phase of matter in the 
inner part of the star will occur. 

Whatever the high density phase, and we assume here it is the quark de- 
confined one, it is softer than the normal phase. Otherwise no phase transition. 
Consequently, a change in the phase of matter will be accompanied by a change 
in the density distribution in the star. In the case of spindown of an isolated 
ms pulsar, the weight of the surrounding part of the star will squeeze the softer 
high density phase that is forming in the core. Conversely, an accreting neutron 
star that is being spun up, will spin out the already present quark phase. In 
either case, the star's moment of inertia will change, and therefore its spin rate 
will change to accommodate the conservation of angular momentum. Timing 
of pulsar and x-ray neutron star rotation is a relatively easy observation, and 
the effect of a phase change in the deep interior of neutron stars on timing or 
frequency distribution is what we study here. 

The deconfinement or reconfinement of quark matter in a rotating star is a 
very slow process because it is governed by the rate of change of period, which is 
very small. This is an advantage. If the processes were fast, we would not likely 
witness the epoch of phase change, which is long, but nevertheless short (but not 
too short) compared to the timescales of spindown of ms pulsars or of spinup 
of accreting neutron stars in low-mass binaries. What we may see in the case of 
isolated ms pulsars is an occasional one that is spinning up, even though losing 
angular momentum to radiation. That would be a spectacular signal. What we 
may see in the case of accreting x-ray neutron stars in low-mass binaries is an 
unusual number of them falling within a small spin-frequency range — the range 
that corresponds to the spinout of the quark matter phase. That also is an easily 
observable signal. For stars of the same mass, the spontaneous spinup of isolated 
pulsars should occur at about the same spin frequency as the clustering in the 
population of neutron star accretors. 

To place the canonical pulsars, ms pulsars, and x-ray accreting neutron stars 
in context, we refer to Fig. Canonical pulsars have large surface magnetic 
fields, of the order of 10^^ to 10^'^ G, and relatively long periods with an average 
of about 0.7 s. The ms pulsars have low fields, of the order of 10^ to 10^ G, 
and periods ranging from 1 to 10 ms. The x-ray stars that are accreting matter 
from a low-mass non-degenerate companion are believed to be the link 
between the populations of canonical and ms pulsars — their path is indicated 
schematically in the figure. (Of course, defining ms pulsars to be those lying in 
the range of periods 1 to 10 ms is arbitrary.) 

At present, fewer ms pulsars have been discovered than canonical pulsars. 
This is likely to be a selection effect. It is much more difficult to detect ms pul- 
sars because there is radio noise in all directions (and therefore a possible signal) 
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Fig. 1. Magnetic field of radio pulsars as a function of rotation period. Death line shows 
a combination of the two properties beyond which the beaming mechanism apparently 
fails. Canonical pulsars are designated as (1), millisecond pulsars as (3) and a schematic 
path for x-ray accreting neutron stars as (2). From Ref. |Q. 



and because of the dispersion by the intersteUar electrons (from ionized hydro- 
gen) of any pulsed signal which might be present in the direction the telescope is 
pointing. If present the signal is weak and contains a band of radio frequencies. 
Because of the different time-delay of each frequency, the frequencies are placed 
into a number of bins. For ms pulsars, the time-delay across the range of frequen- 
cies in a pulse is greater than the interval between pulses. Detection therefore 
depends on a tedious analysis that corrects the time-delay in each frequency 
bin for an assumed density of intervening interstellar gas; this is repeated for a 
succession of assumed densities. If the process converges to a periodic pulse, a 
pulsar has been discovered. Otherwise, the radio telescopes are pointed in a new 
direction. Systematically covering the sky is obviously very costly and time con- 
suming, and moreover, the present instrumentation is unable to detect pulsars 
with periods less than 1 millisecond. 

2 Spontaneous Spinup of Millisecond Pulsar 

We discuss here the effect that a change of phase in the core of a neutron star 
can have on the rotational properties of ms pulsars ||J^,§,|| . The analysis begins 
with the energy-loss equation for a rotating magnetized star whose magnetic 
axis is tilted with respect to the rotation axis. It is of historical interest to note 
that even before pulsars were discovered and soon identified with hypothetical 
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neutron stars, Pacini had postulated the existence of a highly magnetized rotat- 
ing neutron star inside the Crab nebula as the energy source of the nebula and 
inferred some of the star's properties iQ. It was already known before 1967 that 
the nebula, formed in a supernova in 1054, is being accelerated and illuminated 
by a power source amounting to ^ 4 x 10'^* erg/s. This compares with a solar 
luminosity of Lq ~ 4 x 10'^^ erg/s. The power output of the rotating Crab pulsar 
equals that of 100,000 suns. Equating the power input to the nebula with the 
power radiated by a rotating magnetic dipole, 

4 X 10^8 ergs/s - - (^^^') = sin^ a , (1) 

and knowing the period and rate of change of period of the pulsar, 

— s, P - 4 X 10"^^ s/s, (2) 
30 ' ' ^ ' 

the moment of inertia, surface magnetic field strength, and rotational energy of 
the Crab pulsar can be inferred as 

/ - 9 X lO-*"* g cm2 - 70 km^ , 

B - 4 X 10^2 gauss , (3) 
E,ot - 2 X 10^9 ergs - 10^^ MeV . 

For these estimates, we have assumed that sin a = 1, where a is the angle 
between magnetic and rotational axis. (Gravitational units G — c — 1 are used 
frequently. For convenient conversion formula to other units see ch. 3 of Ref. ^.) 

Returning to the effect that a phase change can have on the rotational prop- 
erties of a pulsar through changes in the moment of inertia induced by a phase 
transition, we rewrite Eq. (|^) in greater generality: 

Here we have written for convenience 

C ^{2/?.)R^B^sm^a. (5) 
We find the deceleration equation 

In work previous to ours, /' = dl/df2 was assumed to vanish. This would be a 
good approximation for canonical pulsars but not for millisecond pulsars. 

The angular momentum of a rotating star in General Relativity can be ob- 
tained numerically as a solution of Einstein's equations or as a very complicated 
algebraic expression in a perturbative expansion. The moment of inertia is then 
obtained as / = J / Q. The complication arises in two ways. First, a rotating star 
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Fig. 2. Apparent braking index for pure dipole radiation (n = 3) as function of ro- 
tational frequency for stars of two different masses containing hyperons (but no first 
order phase transition). (The mass indicated is for slow rotation.) 



sets the local inertial frames into rotation. This is referred to as frame dragging. 
Second, the structure of the rotating star depends on the frame dragging fre- 
quency, which is position dependent, and on the spacetime metric, which also 
depends on lo. Algebraic expressions were obtained in Refs. |T]],|9|. 

The index n in equation Eq. equals three for magnetic dipole radiation. 
In principle, it can be measured in terms of the frequency f2 and its first two 
time derivatives. But the dependence of / on i7 introduces a correction so that 
the dimensionless measurable braking index is given by 



niil) ~ —. — 



2i + rn 



(7) 



So the measurable braking index will differ from 3 for ms pulsars, if for no other 
reason than that the centrifugal deformation of the star relaxes as the star spins 
down. Such a dependence on f2 is illustrated in Fig. || 

It is especially noteworthy that the magnitude of the signal — the deviation of 
n{f2) from n — depends on dl / dfl which is large, especially for millisecond pul- 
sars, but the duration of the signal depends on dfi/dt which is small. Therefore, 
the variation of the braking index over time is very slow, the time-scale being 
astronomical. So what one would observe over any observational era is a constant 
braking index n(J?) that is less (if no phase transition) than the n characteristic 
of the energy-loss mechanism Eq. (^) even if the radiation were a pure multipole. 
We will see that the effect of a phase transition on the measurable braking index 
can be (but not necessarily) much more dramatic. (We distinguish between the 
constant n = 3 of the energy-loss equation and the measurable quantity of Eq. 

(0)0 
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We turn now to our original postulate — that canonical pulsars have a quark 
matter phase in their central region but that ms pulsars which are centrifugally 
diluted, do not. In the slow course of spindown the central density will rise how- 
ever. When the critical density is reached, stiff nuclear matter will be replaced 
slowly in the core by highly compressible quark matter. The overlaying layers of 
nuclear matter weigh down on the core and compress it. Its density rises. The 
star shrinks — mass is redistributed with growing concentration at the center. 
The by-now more massive central region gravitationally compresses the outer 
nuclear matter even further, amplifying the effect. The density profile for a star 
at three angular velocities, (1) the limiting Kepler angular velocity at which the 
star is stretched in the equatorial plane and its density is centrally diluted, (2) 
an intermediate angular velocity, and (3) a non-rotating star, are shown in Fig. 
^. We see that the central density rises with decreasing angular velocity by a 
factor of three and the equatorial radius decreases by 30 percent. In contrast, 
for a model for which the phase transition did not take place, the central density 
would change by only a few percent |p^ . The phase boundaries are shown in 
Fig. ^ from the highest rotational frequency to zero rotation. 

The redistribution of mass and shrinkage of the star change its moment of 
inertia and hence the characteristics of its spin behavior. The star must spin up 
to conserve angular momentum which is being carried off only slowly by the weak 
electromagnetic dipole radiation. The star behaves like an ice skater who goes 
into a spin with arms outstretched, is slowly spun down by friction, temporarily 
spins up by pulling the arms inward, after which friction takes over again. 

The behavior of the moment of inertia in the vicinity of the critical region 
of growth of the quark matter core is shown in Fig. ^. The critical period of 
rotation is P = I-k j fl ^ 4.6 ms. Much the same phenomenon was observed in 
rotational nuclei during the 1970s and had been predicted by Mottelson 

and Valatin though of course the phase changes are different (see Fig. ^). 
The mechanism is evidently quite robust, but it need not occur in every model 
star. In the present instance, the mass of the neutron star is very close to the 
limiting mass of the non-rotating counterpart. This does not necessarily mean 
that it will be a rare event. The mass of neutron stars is bounded from below by 
the Chandrasekhar mass limit which is established by electron pressure and has 
a value of \AMq. The limiting mass of neutron stars may be very little more 
because of the softening of the nuclear equation of state by hyperonization and 
quark deconfinemcnt. This could be the reason that neutron star masses seem 
to lie in a very small interval |^ . 

As already pointed out, the progression in time of the growth of the quark 
core is very slow, being governed by the weak processes that cause the loss of 
angular momentum to radiation. Using the computed moment of inertia for a 
star of constant baryon number as shown in Fig. ^ we can integrate Eq. (^) to 
find the epoch of spinup to endure for 2 x 10^ y. This is a small but significant 
fraction of the spin-down time of ms pulsars which is ~ lO^y. So if ms pulsars are 
near their limiting mass and are approximately described by our model, about 
1 in 50 isolated ms pulsars should be spinning up instead of spinning down. 
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Fig. 3. Mass profiles as a function of equatorial radius of a star rotating at three 
different frequencies, as marked. At low frequency the star is very dense in its core, 
having a 4 km central region of highly compressible pure quark matter. At intermediate 
frequency, the pure quark matter phase is absent and the central 8 km is occupied by 
the mixed phase. At higher frequency (nearer fix) the star is relatively dilute in the 
center and ccntrifugally stretched. Inflections at e = 220 and 950 are the boundaries of 
the mixed phase. 

Fig. 4. Radial boundaries at various rotational frequencies separating (1) pure quark 
matter, (2) mixed phase, (3) pure hadronic phase, (4) ionic crust of neutron rich nuclei 
and surface of star. The pure quark phase appears only when the frequency is below 
f2 ~ 1370 rad/s. Note the decreasing radius as the frequency falls. The frequencies of 
two pulsars, the Crab and PSR 1937+21 are marked for reference. 



Presently, about 60 ms pulsars have been identified, and about half of them are 
isolated, the others being in binary systems. Because the period of ms pulsars 
can be measured with an accuracy that rivals atomic clocks, identification of the 
direction of change of period would not take long. For example, PSR1937+21 
has a period (measured on 29 November 1982 at 1903 UT) 

P = 1.5578064487275(3) ms . 

Its rate of change of period is a mere P ~ 10~^^. But because of the high accuracy 
of the period measurement, it would take only two measurements spaced 0.3 
hours apart to detect a unit change in the last significant figure, and hence to 
detect in which direction the period is changing. 

In principal, the period and first two time derivatives can be measured. The 
dimensionless ratio formed from them yields the braking index n of the energy- 
loss mechanism, corrected by a term that depends on changes in the moment 
of inertia. In the case that a phase change causes a spinup of the pulsar, as in 
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Fig. 5. Moment of inertia of a neutron star at angular velocities for which the central 
density rises from below to above critical density for the pure quark matter phase as 
the centrifugal force decreases. Time flows from large to small /. The most arresting 
signal of the phase change is the spontaneous spin-up that an isolated pulsar would 
undergo during the growth in the region of pure quark matter. (Adapted from Ref. 
§■) 

Fig. 6. Nuclear moment of inertia as a function of squared frequency for ^^®Er, show- 
ing backbending in the nuclear case. Quantization of spin yields the unsmooth curve 
compared to the one in Fig. 0. 



Fig. ^, the measurable dimensionless quantity has two singularities at the 

frequencies at which dl /dfl switches between ±00. From Eq. it is clear that 
fi will pass through zero and change sign at both turning points. Therefore, 
the measurable braking index Eq. will have singularities as shown in Fig. 
0. The braking index with a nominal value of 3 will make enormous excursions 
from that value. The braking index is shown in Fig. ^ over one decade in time, 
and has an anomalous value for 10^ y. However, the second time derivative has 
never been measured for a ms pulsar because the rate of change of frequency 
is so slow. So the one, and easily observable signal, is the spontaneous spinup 
of an isolated ms pulsar. Such a pulsar has not yet been observed. But then, 
only about 30 isolated ms pulsars have so far been detected, whereas, according 
to our estimate, only 1 in 50 ms pulsars of mass close to the maximum would 
presently be passing through the epoch in which a quark matter core grows. 

We briefly describe the nuclear and quark matter phases used in this and the 
following sections. More details can be found in the Appendix and in Ref. 
The initial mass of the star in our examples is 1A2Mq. Briefly, confined nuclear 
matter is described by a covariant Lagrangian describing the interaction of the 
members of the baryon octet with scalar, vector and vector-isovector mesons and 
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Fig. 7. Braking index for a pulsar that passes through a change of phase in the central 
region as a function of angular velocity. 

Fig. 8. Braking index as a function of time over one decade in the critical region in 
case of a phase transition in the interior of the star. (From Ref. Q.) 



solved in the meanfield approximation. Quark matter is described by the MIT 
bag model. 

3 Accreting X-ray Neutron Stars in Binaries 

Some neutron stars have non-degenerate companions. Such neutron stars are 
radio silent because a wind from the hot surface of the companion disperses 
the pulsed radio signal which a rotating magnetized neutron star would otherise 
radiate into space. Late in the life of the neutron star, when the slowly evolving 
companion begins to overflow the Roche lobe, mass transfer onto the neutron 
star commences. The drag of the magnetic dipole torque will be eclipsed by the 
transfer of mass and angular momentum onto the neutron star. It has begun its 
evolution from an old slowly rotating neutron star with long period and high 
magnetic field, to a ms pulsar with low field (sometimes referred to as a "recycled 
pulsar"). During the long intermediate stage, when the surface and accretion ring 
are heated to high temperature, the star radiates x-rays. 

Any asymmetry in the mass accretion pattern will cause a variability in x-ray 
emission. At an accretion rate of lO^^Afo/y it would take only 10^ y to spin 
up the neutron star to a period of 2 ms (500 Hz). Consequently, millisecond 
variability in the x-ray luminosity is expected and observed. 

Accreting x-ray neutron stars provide a very interesting contrast to the spin- 
down of isolated ms pulsars [|6|. Presumably they are the link between the 
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canonical pulsars with mean period of 0.7 sec and the ms pulsars If 
the critical deconfinement density falls within the range spanned by canonical 
pulsars, quark matter will already exist in them but may be "spun" out of x-ray 
stars as their rotational frequency increases during accretion. We can anticipate 
that in a certain frequency range, the changing radial extent of the quark mat- 
ter phase will actually inhibit changes in frequency because of the increase in 
moment of inertia occasioned by the gradual disappearance of the quark matter 
phase. Accreters will tend to spend a greater length of time in the critical fre- 
quencies than otherwise. There will be an anomalous number of accreters that 
appear at or near the same frequency. This is what was found recently in data 
obtained with the Rossi X-ray Timing Explorer (RXTE). For an extensive re- 
view of the discoveries made in the short time since this satellite was launched 
(1995), see Ref. |l|]. 

The spinup evolution of an accreting neutron star is a more complicated 
problem than that of the spindown of an isolated ms pulsar of constant baryon 
number. It is complicated by the accretion of matter (M ^ 10^ ^^M© yr~^), a 
changing magnetic field strength (from B ~ 10^^ to ~ 10^ G), and the interac- 
tion of the field with the accretion disk. 

The change in moment of inertia as a function of rotational frequency caused 
by spinup due to accretion is similar to that described in the previous section, 
but in reverse However, there are additional phenomena as just mentioned. 
The spin-up torque of the accreting matter causes a change in the star's angular 
momentum according to the relation [ p"8Hl9| , |20| 

^J=^(W)^NA{rrn)-NMirc). (8) 

at at 

The first term on the right-hand-side is the torque exerted on the star by a 
mass element rotating at the base of the accretion disk with Keplerian velocity 
lok- Denoting this distance with r,,!, one readily finds that -/VA(^m) is given by 
(G = c=l) 

NA{rm) = rj^n MUJK 

= M = M , (9) 

where M stands for the accretion rate, and l^r^^) is the specific angular momen- 
tum of the accreting matter (angular momentum added to the star per unit mass 
of accreted matter). The second term on the right-hand-side of Eq. (^) stands 
for the magnetic plus viscous torque term (k ~ 0.1), 

NMirc)=Kfi^r-^, (10) 

with fi = R^B the star's magnetic moment. Upon substituting Eqs. (||) and (^0|) 
into (^ and writing the time derivative d{If2)/dt as d{IQ)/dt ~ {dl/dt)[2 + 
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I{dQ /dt), the time evolution equation for the angular velocity fl of the accreting 
star can be written as 

m'^ = mt) -^m' rsr' ■ (n) 

The quantities rm and Tc denote fundamental length scales of the system. The 
former, as mentioned above, denotes the radius of the inner edge of the accretion 
disk and is given by ~ 1) 

rm^CrA. (12) 
The latter stands for the co-rotating radius defined as 

r,^{MQ~Y'- (13) 

Accretion will be inhibited by a centrifugal barrier if the neutron star's magne- 
tosphere rotates faster than the Kepler frequency at the magnetosphere. Hence 
''m < ''c, otherwise accretion onto the star will cease. A further fundamental 
lengthscale is set by the Alfen radius ta, which enters in Eq. It is the 

radius at which the magnetic energy density, i?^(r)/87r, equals the total kinetic 
energy density, p(r)v'^(r)/2, of the accreting matter. For a dipole magnetic field 
outside the star of magnitude 

Bir) = ^ , (14) 

this condition reads 

^S2(rA) = ip(rA)«'(rA), (15) 

where w(rA) is of the order of the Keplerian velocity, or, which is similar, the 
escape velocity at distance ta from the neutron star, 



virj,)^(^y\ (16) 



The density p{rA) in Eq. (|T^) can be replaced by 

, ^ M 

^ ''A = -j — — r , (17) 

which follows from the equation of continuity. With the aid of Eqs. (|l^), ( p^ 
and (p^), one obtains from (15) for the Alfen radius the following expression: 

^«=f^)"'^ (18) 

V2MM2/ 



It is instructive to write this equation as 

13 ,V-2/7 ,4/7 



^1/7 



r^ = 7xlO'Aa^' pZ'[^] km, (19) 
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where M_io = M/(10"i° M© yr^^) and /iso = /i/(10^° G cm^). It shows that 
canonical accreters (M_io — 1) with strong magnetic fields of B ~ 10^^ G have 
Alfen radii, and thus accretion disks, that are thousands of kilometers away from 
their surfaces. This is dramatically different for accreters whose magnetic fields 
have weakened over time to values of ~ 10* G, for instance. In this case, the 
Alfen radius has shrunk to ~ 40 km, which is just a few times the stellar radius. 
We assume that the magnetic field evolves according to 

= B(oo) + [5(0) - B(oo)]e-*/*<i (20) 

with t = at the start of accretion, and where i?(0) = 10^^ G, B{oo) = 10* G, 
and = 10^ yr. Such a decay to an asymptotic value seems to be a feature of 
some treatments of the magnetic field evolution of accreting neutron stars |^ . 
Moreover, it expresses the fact that canonical neutron stars have high magnetic 
fields and ms pulsars have low fields (see Fig. |l|). Beyond that, as just mentioned, 
the condition that accretion can occur demands that < which inequality 
places an upper limit on the magnitude of the magnetic field of ms neutron star 
accreters of 2 to 6 x 10* G @. 

Frequently, it has been assumed that the moment of inertia in Eq. (|l^) does 
not respond to changes in the centrifugal force, and in that case, the above 
equation yields a well-known estimate of the period to which a star can be 
spun up 1^. The approximation is true for slow rotation. However, the response 
of the star to rotation becomes increasingly important as the star is spun up. 
Not only do changes in the distribution of matter occur but internal changes 
in composition occur also because of changes induced in the central density by 
centrifugal dilution |^ ; both changes effect the moment of inertia and hence the 
response of the star to accretion. 

The moment of inertia of ms pulsars or of neutron star accreters has to be 
computed in GR without making the usual assumption of slow rotation p^p^ . 
We use a previously obtained expression for the moment of inertia of a rotat- 
ing star, good to second order in Q [ pT| . The expression is too cumbersome to 
reproduce here. Stars that are spun up to high frequencies close to the breakup 
limit (Kepler frequency) undergo dramatic interior changes; the central density 
may change by a factor of four or so over that of a slowly rotating star if a phase 
change occurs during spin- up (cf. Fig. ||) |0,|2^ . 

Figure |^ shows how the moment of inertia changes for neutron stars in binary 
systems that are spun up by mass accretion according to Eq. (|l|) until QAMq 
has been accreted. The neutron star models are fully described in Ref. ^ and 
references therein and briefly in the Appendix of this article. The initial mass 
of the star in our examples is 1.42M0. In one case, it is assumed that a phase 
transition between quark matter and confined hadronic matter occurs, and in 
the other that it does not. This accounts for the different initial moments of 
inertia, and also, as we see, the response to spinup. Three accretion rates are 
assumed, which range from M_io = 1 to 100 (where M_io is measured in units 
of 10~^°Mq per year). These rates are in accord with observations made on low- 
mass X-ray binaries (LMXBs) observed with the Rossi X-ray Timing Explorer 
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Fig. 9. Moment of inertia of neutron stars as a function of time witii (solid curves) 
and witiiout (dashed curves) quark matter core assuming 0.4Mq is accreted. Results 
for three average accretion rates are illustrated. 



[ p^ . The observed objects, which are divided into Z sources and A(toll) sources, 
appear to accrete at rates of M-iq ^ 200 and Af_io ~ 2, respectively. Although 
in a given binary, M varies on a timescale of days, we take it to be the constant 
average rate in our calculations. 

Figure |l^ shows the spin evolution of accreting neutron stars as determined 
by the changing moment of inertia and the spin evolution equation, Eq. (pi]). 
Neutron stars witiiout quark matter in their centers are spun up along the dashed 
lines to equilibrium frequencies between about 600 Hz and 850 Hz, depending on 
accretion rate and magnetic field. The dl/dt term for these sequences manifests 
itself only insofar as it limits the equilibrium periods to values smaller than the 
Kepler frequency, vk- In both Figs. ^ and |l^ we assume that OAMq is accreted. 
Otherwise, the maximum frequency attained is less, the less matter is accreted. 

The spin-up scenario is dramatically different for neutron stars in which a first 
order phase transition occurs. In this known from Fig. ^, the temporal 

conversion of quark matter into its mixed phase of quarks and confined hadrons 
is accompanied by a pronounced increase of the stellar moment of inertia. This 
increase contributes so significantly to the torque term N{rc) in Eq. (^l]) that the 
spinup rate dfl / dt is driven to a plateau around those frequencies at which the 
pure quark matter core in the center of the neutron star gives way to the mixed 
phase of confined hadronic matter and quark matter. The star resumes ordinary 
spin-up when this transition is completed. The epoch during which the spin rates 
reach a plateau are determined by attributes like the accretion rate, magnetic 
field, and its assumed decay time. The epoch lasts between ^ 10^ and 10^ yr 
depending on the accretion rate at the values taken for the other factors. 
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Fig. 10. Evolution of spin frequencies = il/2n) of accreting neutron stars with 
(solid curves) and without (dashed curves) quark deconfinement if 0.4Mq is accreted. 
(If the mass of the donor star is less, then so is the maximum attainable frequency.) 
The spin plateau around 200 Hz signals the ongoing process of quark confinement in 
the stellar centers. Note that an equilibrium spin is eventually reached which is less 
than the Kepler frequency. 



We can translate the information in Fig. |l^ into a frequency distribution 
of X-ray stars by assuming that neutron stars begin their accretion evolution 
at the average rate of one per million years. A different rate will only shift 
some neutron stars from one bin to an adjacent one, but will not change the 
basic form of the distribution. The donor masses in the binaries are believed to 
range between 0.1 and O.4Af0. For lack of more precise knowledge, we assume a 
uniform distribution of donor masses (or mass accreted) in this range and repeat 
the calculation shown in Fig. 10 at intervals of O.lAf©. 



The result for the computed distribution of rotational frequency of x-ray neu- 
tron stars is shown in Fig. it is striking. Spinout of the quark matter core 
as the neutron star spins up is signalled by a spike in the distribution which 
would be absent if there were no phase transition in our model of the neutron 
star. We stress that what we plot is our prediction of the relative frequency dis- 
tribution of the underlying population of x-ray neutron stars — but the weight 
given to the spike as compared to the high frequency tail depends sensitively 
on the weight with which the donor masses are assigned and the initial mass 
distribution of neutron stars in LMXBs. As already mentioned, we give equal 
weight to donor masses between 0.1 and O.4M0 and the initial mass of the neu- 
tron star is 1.42M0. The objects above 400 Hz in Fig. |ri| are actually unstable 
and will collapse to black holes. Donors of all masses in the range just mentioned 
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Fig. 11. Frequency distribution of X-ray neutron stars. Calculated distribution (open 
histogram) for the underlying population is normalized to the number of observed 
objects (18) at the peak. (The normalization causes a fractional number to appear in 
many bins of the calculated distribution.) Data on neutron stars in LMXBs (shaded 
histogram) is from Ref. [^]. See text and especially the reference for caveats to the 
interpretation. The spike in the calculated distribution corresponds to the spinout 
of the quark matter phase and the corresponding growth of the moment of inertia 
as compressible quark matter is replaced by relatively incompressible nuclear matter. 
Otherwise the spike would be absent. 



contribute to neutron stars of spin up to 400 Hz. Neutron stars of lower initial 
mass than our 1.42Mq with donors of mass at the higher end of their range 
will produce spins above 400 Hz. In other words, the relative population in the 
peak as compared to the background will be sensitive to the unknown factors (1) 
accretion rate (2) initial mass distribution of neutron stars in LMXBs (3) mass 
distribution of donor stars. However, the position of the peak in the spin distri- 
bution of x-ray neutron stars is a property of nuclear matter and independent 
of the above unknowns. 

The calculated concentration in frequency of x-ray neutron stars is centered 
around 200 Hz; this is about 100 Hz lower than the observed spinup anomaly 
(see discussion below). This discrepancy should not be surprising in view of our 
total ignorance]^ of the equation of state above saturation density of nuclear 
matter and the necessarily crude representation of hadronic matter in the two 
phases in the absence of relevant solutions to the fundamental QCD theory of 
strong interactions. We represent the confined phase by relativistic nuclear field 
theory and the deconfined phase by the MIT bag model. However crude these or 
any other models of hadronic matter may be, the physics underlying the effect 
of a phase transition on spin rate is robust, although not inevitable. We have 

^ There are upcoming radioactive beam experiments from which it is hoped to gain 
information on the equation of state of asymmetric nuclear matter 
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Fig. 12. Data on frequency distribution of millisecond pulsars (1 < P < 10 ms). 
Frequency bins are 50 Hz wide. 



cited the example of an analogous phenomenon found in rotating nuclei in the 
previous section. 

The data that we have plotted in Fig. ^is gathered from Tables 2, 3, and 4 
of the review article of van der Klis concerning discoveries made with the Rossi 
X-ray Timing Explorer, launched near the end of 1995 [0. The interpretation 
of millisecond oscillations in the x-ray emission, either that found in bursts or of 
the difference between two observed quasi-periodic oscillations (QPOs) in x-ray 
brightness, is ambiguous in some cases. In particular, the highest frequency near 
600 Hz in the "observed" data displayed in our figure, may actually be twice 
the rotational frequency of the star p^ . The millisecond variability in x-ray 
phenomena associated with accretion onto neutron stars that has been observed 
since the launch of the satellite was anticipated several decades ago. However, 
the consistency of the phenomena from one binary to another raises questions 
about interpretation. In this sense, the field is quite young. We refer to the above 
cited review article for details and references to the extensive literature. 

Nevertheless, the basic feature will probably survive — a clustering in rota- 
tional frequencies of x-ray neutron stars and a higher frequency tail. Certainly 
there are high frequency pulsars. A histogram of ms pulsar frequencies shows a 
broad distribution around 200 Hz, and a tail extending to ^ 600 Hz as shown 
in Fig. |l^. So both the (sparse) data on X-ray objects and on ms pulsars seem 
to agree on a peak in the number of stars at moderately high rotational fre- 
quency and on an attenuation at high frequency. For ms pulsars, however, the 
attenuation at high frequency may be partly a selection effect due to interstellar 
disperion of the radio signal. 

There have been other suggestions as to the cause of the spike, several of 
which we cite (c.f. [ p8|j29| , |30t ) . These works are concerned with the balance of 
the spinup torque by a gravitational radiation torque. Our proposal has several 
merits: (1) the mechanism involving a change in moment of inertia triggered by 
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a phase transition in the stellar core due to changing density profile in the star 
as its spin changes is robust — it is known to occur in rotating nuclei; (2) The 
phase transition should occur in reverse in isolated ms pulsars and in neutron 
star accretors and at about the same frequency for similar mass stars; (3) The 
phase transition causes accretion induced spinup to stall for a long epoch, but 
to resume after the quark core has been expelled, thus accounting for high spin 
objects like very fast ms pulsars as well as the clustering in spin of the popuation 
of accretors. 

4 Evolution from Canonical to Millisecond Pulsar 

In the foregoing we have discussed possible signals of a phase transition in iso- 
lated ms pulsars and in accreting x-ray neutron stars in binary orbit with a 
low-mass companion. Spinup by mass accretion is believed to be the pathway 
from the relatively slowly rotating canonical pulsars formed by conservation of 
angular momentum in the core collapse of massive stars, and the rapidly rotating 
millisecond pulsars jljj^ . In this section we trace some of the possible evolution- 
ary routes under the various physical conditions under which accretion occurs 

[ill- . 

The evolutionary track between canonical and ms pulsars in the coordinates 
of magnetic field strength and rotational period (refer to Fig. ^ will depend on 
the rate of mass accretion, its duration, the centrifugal change in the moment 
of inertia of the star, the strength of the magnetic field and timescale of its 
decay, and possibly other affects. Many papers have been devoted to the decay 
of the magnetic field. It is an extremely complicated subject, with many physical 
uncertainties such as the actual location of the field, whether in the core or crust, 
the degree to which the crust is impregnated with impurities, crustal heating and 
resultant reduction in conductivity and therefore increase in ohmic field decay, 
screening of the magnetic field by accreted material, and so on. 

The field is believed to decay only weakly due to ohmic resistance in canonical 
pulsars, but very significantly if in binary orbit with a low-mass non-degenerate 
star, when the companion fills its Roche lobe. This era can last up to 10^ y and 
cause field decay by several orders of magnitude. For a review of the literature 
and several evolutionary scenarios, see Refs. p0| , ^l]j3^ , |3^ , ^ . 

While there is no consensus concerning the magnetic field decay, observa- 
tionally, we know that canonical pulsars have fields of ~ 10^^ to lO^^G, while 
millisecond pulsars have fields that lie in the range ~ 10^ to lO^G. We shall rely 
on this observational fact, and assume that the field decays according to Eq. 
(|0|). where B{oo) = 10^ G, B{0) = lO^^ G and U = 10^ to 10^ yr. Moreover, 
this is the general form found in some scenarios [^. However, we shall also make 
a comparison with a purely exponential decay. 

There are three distinct aspects to developing an evolutionary framework. 
One has to do with the accretion process itself, which has been developed by a 
number of authors in the framework of classical physics (||l^,|l9|;^) and which we 
employed in the previous section. Another has to do with the field decay, which 
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in a complete theory will be coupled to the accretion process. The third aspect 
has to do with the structure of the neutron star and its response to added mass, 
but most especially to its response to changes in rotational frequency due to the 
changing centrifugal forces. 

Typically, the moment of inertia has been computed in general relativity 
for a non-rotating star |22 2^. It is based on the Oppenheimer-Volkoff metric. 
However, for the purpose of tracing the evolution of an accreting star from ~ 1 
Hz to 400 — 600 Hz we do not neglect the response of the shape, structure 
and composition of the star as it is spun up over this vast range of frequencies 
from essentially zero to values that approach the Kepler frequency. Nor do we 
neglect the dragging of local inertial frames. These features are included in our 
calculation of the tracks of neutron stars from canonical objects starting with 
large fields and very small frequencies at the "death line" to the small fields but 
rapid rotation of millisecond pulsars. However, the expression of the moment of 
inertia and a definition of the various factors that enter are too long to reproduce 
here. We refer instead to our derivation given in Refs. which is computed 

to second order in the rotational angular momentum and found to be accurate 
to ^ 10 % when compared to numerical solutions for a rotating star [p5| . In the 
present context, numerical solutions were obtained in Ref. [Bm, and semianalytic 



approximations were employed in Ref. |37 




Fig. 13. Evolutionary tracks traced by neutron stars in the X-ray accretion stage, 
beginning on the death line with large B field and ending as millisecond stars, for 
various accretion rates. Here = lO^yr. 
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The initial conditions for the evolution are arbitrary to a high degree. Canon- 
ical pulsars have a broad range of magnetic field strengths. The period of the 
pulsar at the time that the companion overflows its Roche lobe and accretion 
commences is also arbitrary. Any observed sample of x-ray accreters presumably 
spans a range in these variables. For concreteness, we assume that the pulsar 
has a field of 10^^ G and that the period of the pulsar is 1 ms when accretion 
begins. The donor mass in the low-mass binaries are in the range 0.1 to OAMq. 
Our sample calculations are for accretion of up to OAMq. 
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Fig. 14. Evolutionary tracks traced by neutron stars in the X-ray accretion stage, 
beginning on the death line with purely exponential decay of the B field. As in Fig. 
= lOVr. 



We find essentially a continuum of evolutionary tracks in the B ~ P plane 
according to the rate at which matter is accreted from the companion, and 
the rate at which the magnetic field decays. The evolutionary tracks essentially 
fill all the space in the B-P plane, starting at our assumed initial condition 
of an old canonical pulsar with field of 10^^ G, and extending downward in 
field strength, broadening to fill the space on both sides of the deathline, and 
extending to the small periods of millisecond pulsars. All are potential tracks 
of some particular binary pair, since accretion rates vary by several orders of 
magnitude and presumably so do decay rates of the magnetic field. As a first 
orientation as to our results and how they relate to known pulsars as regards their 
magnetic field strength and their rotational period, we show the evolutionary 
tracks for four different accretion rates given in units of M = 10~^° solar masses 
per year in Fig. Ol The decay rate of the magnetic field is taken to have the 
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value id — 10® yr in each case. The x-ray neutron star gains angular momentum 
and its period decreases, and over a longer timescale, the magnetic field decays. 
One can see already that a wide swathe of B and P is traced out. 

In the above example, the field was assumed to decay to a finite asymp- 
totic value of 10^ Gauss. A very different assumption, namely that the field 
decays eventually to zero, B{t) = i3(0)e~*/*'i, modifies only the results below 
the asymptotic value, as is seen by comparing Fig. |l^ and However, the con- 
clusion concerning the origin of millisecond pulsars is quite different. For purely 
exponential decay, one would conclude that high frequency pulsars are created 
only in high accretion rate binaries. 

In the remainder of the paper, we assume the field decays to an asymptotic 
value, since from the above comparison we see how exponential decay would 
modify the picture. 

We show time tags on a sample track in Fig. |l^ which provides some sense of 
time lapse. The first part of a track is traversed in short time, but the remainder 
ever more slowly. This shows up also in dP/dt as a function of time. For each of 




Fig. 15. Time tags expressed in years are shown for one of the evolutionary tracks 
corresponding to a decay constant for the magnetic field of td ~ 10^ y. 

Fig. 16. Evolutionary tracks for a neutron star starting at the death line and evolving 
by accretion to lower field and high frequency for an accretion rate 1 in units of 10~^^ 
solar masses per year, and for three values of the magnetic field decay rate td as marked. 



three accretion rates we show the dependence on three field decay constants in 
Figs. |l^ and 18. Depending on decay rate of the field and accretion rate, an 
X-ray neutron star may spend some time on either side of the death line, but if 
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it accreted long enough, always ends up as a candidate for a millisecond pulsar if 
the magnetic filed decays to an asymptotic value such as was assumed. However, 
if the field decays exponentially to zero, only high accretion rates would lead to 
millisecond pulsars. Of course, if accretion turns off at some time, the evolution 
is arrested. 




Fig. 17. Evolutionary tracks similar to Fig. y_6| but with a different accretion rate 
{dM/dt)-io = 10. 

Fig. 18. Evolutionary tracks similar to Fig. ^ but with a different accretion rate 
{dM/dt)-w = 100. 



In summary, we have computed the evolutionary tracks in the B — P plane 
due to mass accretion onto neutron stars beginning at the death line with a 
typical field strength of 10^^ Gauss, to shorter periods and low fields. According 
to the assumed accretion rate and decay constant for the magnetic field, the 
tracks indicate that the individual binaries with characteristics ranging from 
Z to Atoll sources will evolve along paths that cover a broad swathe in the 
B — P plane. These include tracks of X-ray stars corresponding to low accretion 
rates that follow a path beyond the death line in the so-called "graveyard'. 
We have assumed two particular forms for the law of decay of the magnetic 
field. (1) The field approaches an asymptotic value of 10® Gauss such as is 
typical of millisecond pulsars. This assumption leads to a particular form for the 
termination of evolutionary tracks. All accreters, no matter what the accretion 
rate, will end with millisecond periods, unless accretion ceases beforehand. (2) If 
instead, we had assumed a purely exponential decay, the tracks would not tend 
to an asymptotic value, but would continue to decrease in the strength of B. 
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The tracks would still cover a broad swathe in the B — P plane. But one would 
conclude that only the higher accretion rate binaries, particularly the Z-sources, 
could produce millisecond period neutron stars. If accretion continues for too 
long a time, the neutron star will be carried to very low fields and across the 
death line, or an overcritical mass will have been accreted, leading instead to a 
black hole. 

5 Appendices 

5.1 First Order Phase Transition in Neutron Stars 

We briefly recall some of the main characteristics of a first order phase transition 
in any substance having more than one conserved quantity, or as we will call it 
in the context of physics — conserved charge — such as baryon number or electric 
charge [Q. These are the two conserved quantities relevant to neutron stars. 
They have zero net electric charge and are made from baryons. 

What makes a substance having more than one conserved charge different 
from a substance having only one, is that charges in the first case can be ex- 
changed between two phases of the substance in equilibrium so as to minimize 
the energy. And the concentration of the charges can readjust at each propor- 
tion of the phases to minimize the energy. There is no such degree of freedom 
in a single-component substance, and there are n — 1 degrees of freedom in a 
substance having n independent components or conserved charges. 

In a single-component substance, the pressure at constant temperature re- 
mains unaltered for all proportions of the two phases as the substance is com- 
pressed. Only the proportion of the phases changes. This is not so for a substance 
of more than one conserved charge. The degree(s) of freedom to readjust the con- 
centrations of the charges at each proportion of the phases causes a change in 
all internal properties of the two phases as their proportion changes including 
their common pressure under conditions of constant temperature. 

The fact that the pressure varies as the proportion of phases in equilibrium 
in neutron star matter has an immediate consequence. The mixed phase of two 
phases in equilibrium will span a finite radial distance in a star. If the pressure 
were a constant for all proportions, then the mixed phase would be squeezed 
out of the star because the pressure varies monotonically, being greatest at the 
center and zero at the edge of the star. In early work on phase transitions in 
neutron stars, they were always treated as having constant pressure at the zero 
temperature of the star, so the mixed phase never appeared in any of those 
models. 

The difference in properties of first order phase transitions in a one-component 
substance and one with several independent components or conserved charges is 
easily proven. First consider a one-component substance for which equilibrium 
of two phases, A and B is expressed by 
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At constant T, the solution for the chemical potential is obviously unique and 
independent of the proportion of the phases. So all properties of the two phases 
remain unaltered as long as they are in equilibrium, no matter the proportion. 

Now consider a substance having two conserved charges (or independent com- 
ponents). For definiteness we consider the two independent conserved charges of 
neutron star matter, baryon number and electric charge, whose densities we de- 
note by p and q. For the corresponding chemical potentials we choose those of 
the neutron, /i„, and electron, /ig. The chemical potentials of all other particles 
can be written in terms of these independent ones. Gibbs phase equilibrium be- 
tween the confined hadronic phase C and the deconfined quark matter phase D 
is now expressed as 

Pc{f^n,f^e,T)=PD{f^,i,fie,T). (22) 

This equation is insufficient to find the chemical potentials. At fixed T, it must 
be supplemented by another, say a statement of the conservation of one of the 
conserved charges. How should that statement be made? If one of the charges 
is the electric charge, demanding that the electric charge density should vanish 
identically in both phases would satisfy the condition of charge neutrality as 
required of a star.Q That is in fact how charge neutrality in neutron stars was 
enforced for many years. However, it is overly restrictive. The net charge must 
vanish, but the charge density need not. Only the integrated electric charge 
density must vanish, / q{r)(Pr = 0.^ We refer to this as global conservation 
rather than local. It releases a degree of freedom that the physical system can 
exploit to find the minimum energy. To express this explicitly, we note that 
according to the preparation of the system, whether in the laboratory, or in a 
supernova, the concentration of the charges is fixed when the system is in a 
single phase. Denote the concentration by 

c = Q/B . (23) 

However, when conditions of temperature or pressure change to bring the system 
into two phases equilibrium, the concentrations in each can be different 

cc = Qc/Bc, CD = Qd/Bd (24) 

provided only that the total charges are conserved, 

Qc + Qd=Q, Bc + Bd=B. (25) 

Here Q and B denote the total electric and baryon charge in a volume V. The 
rearrangement of charges will take place to minimize the energy of the system. 
The force that is responsible for exploiting this degree of freedom in neutron 

^ The Coulomb force is so much stronger than the gravitational that the net charge 

per baryon has to be less than 10"'^^ which we can call zero. 
^ Familiar examples of neutral systems that have finite charge densities of opposite 

sign are atoms and neutrons. 
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stars is the one responsible for the symmetry energy in nuclei and the valley of 
beta stability. Since neutron stars are far from symmetry, the symmetry energy 
is quite large; the difference in Fermi energies of neutron and proton, and the 
coupling of baryon isospin to the neutral rho meson are responsible. 

For a uniform medium, and every sufficiently small region y in a neutron 
star is uniform to high accuracy, the statement of global conservation is 

/ qc{r)d\+ I qD{r)d^r^Vcqc + VDqD = Q (26) 
Jvc Jvo 

where Vc and Vd denote the volume occupied by the two phases respectively. 
This can be written more conveniently as 

(1 - x)'7c(Ain,Me) +Xto(MTi,Aie) ^Q/V = 1 (27) 

where qc denotes the density of the conserved charge in the confined phase, qo 
in the deconfined phase, and 

X = Vd/V, V^Vc + Vd (28) 

is the volume proportion of phase D and q is the volume averaged electric charge 
(which for a star is zero). Now, equations (^) and (|2^) arc sufficient to find 
and /ig. But notice that the solutions depend on the volume proportion x- 
Therefore, also all properties of the two phases depend on their proportion, 
including the common pressure. Having solved for the chemical potentials (and 
all field quantities specified by their equations of motion), the densities of the 
baryon conserved charge in the phases C and D, are given by pc{fJ-m fJ'e) and 
PoifJ-m l^-e)- The volume average of the baryon density is given by an equation 
corresponding to Eq. (p7|). 

Let us now discuss the consequences of opening the degree of freedom embod- 
ied in Eq. (p7|), ie., in allowing electric charge (and strangeness) to be exchanged 
by the two phases in equilibrium so as to achieve the minimum energy at the 
corresponding baryon density. Because of the long range of the Coulomb force 
the Coulomb energy will be minimized when regions of like charge are small, 
whereas the surface interface energy will be minimized when the surface areas 
of the regions of the two phases is small. These are opposing tendencies, and in 
first order can be reconciled by minimizing their sum. A Coulomb lattice will 
of such a size and geometry of the rare phase immersed at spacings 
1C| in the dominant phase so as to minimize the energy. 
In better approximation, the total energy, consisting of the sum of bulk en- 
ergies, the surface and Coulomb energies and higher corrections such as the 
curvature energy, 

^^Total ~ -^Bulk + -Esurfaco + -^Coulomb + -^Curvature + ' ' ' , (29) 

should be minimized. In still better approximation, the convenient partition of 
the energy as above would be replaced by a lattice calculation of the total energy 
[El]. In general, the opening of the degree(s) of freedom to conserve charges 
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Fig. 19. Equation of state for a first order pliase transition from neutron star matter in 
its confined to deconfined phases (marked hybrid). Note the monotonically increasing 
pressure. Comparison is made with the equation of state of neutron star matter in the 
confined phase with nucleons hyperons and leptons in equilibrium. 

globally rather than locally, or any other arbitrary way, can only lower the total 
energy from the value it would have were the degree of freedom closed, or in 
very special cases leave it unchanged. However, in the case of a neutron star, 
the degree of freedom allows the bulk energy in the normal phase to be lowered 
by decreasing the charge asymmetry of neutron star matter, so exchange of 
charge between the confined and deconfined phases is evidently favorable. It is 
unphysical to choose an arbitrary value of surface tension such that the mixed 
phase is energetically unfavored. The surface tension should be calculated self- 
consistently by minimizing the total energy Eq. (|2^) , when possible . 

However, our purpose here is not to calculate the geometric structure^ of the 
mixed phase but to exhibit the equation of state for a first order deconfinement 
transition in neutron star matter according to the above principles. Fig. |l9| shows 
the equation of state with the pure phases at low and high density and the mixed 
phase between. Note, as pointed out before that though the phase transition is 
first order, having as it does a mixed phase, the pressure (and all other internal 
properties) vary with density. The monotonic increase of pressure contrasts with 
early treatments of phase transitions in neutron stars prior to 1991, for which 
approximations rendered a constant pressure in the mixed phase [^^. For a 
constant pressure phase transition, gravity will squeeze out the mixed phase. 
Pressure is a monotonically decreasing function of distance from the center of a 
star, just as it is in our atmosphere. 



See Ref. |47]] for a review. 
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5.2 Description of the Confined and Deconfined Phases 

To describe the confined pliase of neutron star matter, we use a generalization 
of relativistic nuclear field theory solved at the mean field level in which 
nucleons and hyperons (the baryon octet) are coupled to scalar, vector and 
vector-isovector mesons. A full description of how the theory can be solved for 
neutron star matter can be found in PJ44| . The Lagrangian is 

B 

-^9pBlt.T ■ p^')'ipB + ^{df^ad^a - ml<j^) 

+ E ^x{nf^9''-mx)^x. (30) 

The summ on B is over all charege states of the baryon octet. The parameters of 
the nuclear Lagrangian can be algebaically determined (see 2'nd ed. of Q) so that 
symmetric nuclear matter has the following properties: binding energy of sym- 
metric nuclear matter B/A = —16.3 MeV, saturation density p — 0.153 fm^'^, 
compression modulus K — 300 MeV, symmetry energy coefficient Osym = 32.5 
MeV, nucleon effective mass at saturation m*^^ = 0.7m and ratio of hyperon 
to nucleon couplings = 0.6, x^j = 0.653 = Xp that yield, together with the 
foregoing parameters, the correct A binding in nuclear matter ]4^ ). 

Quark matter is treated in a version of the MIT bag model with the three 
light flavor quarks (m„ = rn^ = 0, = 150 MeV) as described in Ref. Q. A 
value of the bag constant i?^/^ = 180 MeV is employed. 
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